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Coherence and Rydberg blockade of atomic ensemble qubits 
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We demonstrate \W) state encoding of multi-atom ensemble qubits. Using optically trapped Rb 
atoms the T 2 coherence time is 2.6(3) ms for N = 7.6 atoms and scales approximately inversely with 
the number of atoms. Strong Rydberg blockade between two ensemble qubits is demonstrated with 
a hdelity of 0.89(1) and a hdelity of ~1.0 when postselected on control ensemble excitation. These 
results are a significant step towards deterministic entanglement of atomic ensembles. 

PACS numbers: 03.67.-a, 42.50.Dv, 32.80.Rm 


Qubits encoded in hyperfine states of neutral atoms 
are a promising approach for scalable implementation of 
quantum information processing[lj. While a qubit can be 
encoded in a pair of ground states of a single atom, it is 
also possible to encode a qubit, or even multiple qubits, 
in an N atom ensemble by using Rydberg blockade to 
enforce single excitation of one of the qubit states^, Q. 
Ensemble qubits have several interesting features in com¬ 
parison to single atom qubits. Using an array of traps it 
is simpler to prepare many ensemble qubits with TV > 1 
for each ensemble, than it is to prepare an array with 
exactly one atom in each trap which remains an out¬ 
standing challenge 0-0. In addition, a \W) state en¬ 
semble qubit encoding is maximally robust against loss 
of a single atomQ, which can be remedied with er¬ 
ror correction protocols while atom loss is a critical 
error for single atom qubits. Furthermore an ensem¬ 
ble encoding facilitates strong coupling between atoms 
and light, an essential ingredient for quantum network¬ 
ing protocols @ and atomic control of photonic interac¬ 
tions in Rydberg blockaded ensembles [l0|. As the atom- 
light coupling strength grows with the number of atoms, 
recent exDeriments[l(ij|. 111 and theory proposals [l^ are 
based on ensembles with N > 100. We are focused here 
on studying the physics of ensembles for computational 
qubits and therefore work with smaller ensembles with 
up to ~ 10 atoms. 


In this letter we demonstrate and study the coherence 
and interactions of atomic ensemble qubits. We measure 
the T 2 coherence time of ensemble qubits achieving a ra¬ 
tio of coherence time to single qubit tt rotation time of 
^ 2600. We furthermore proceed to demonstrate strong 
Rydberg blockade between two, spatially separated en¬ 
semble qubits. Together with the recent demonstration 
of entanglement between a Rydberg excited ensemble and 
a propagating photon[l^ these results establish a path 
towards both local and remote entanglement of arrays 
of ensemble qubits, which will enable enhanced quantum 
repeater architectures [l^ . 


The computational basis states of the ensemble qubits 


are 

1 ^ 

|0) = |0i...0^), |I) = ^^|0i02...1,...0^), (1) 

where |0j) and |lj) are two ground states of the atom 
in an N atom sample [l5j|. The state |I), which is a sym¬ 
metric superposition of one of the N atoms being excited, 
is commonly referred to as a \W) state in the quantum 
information literature. 

Gate protocols for ensemble qubits differ slightly from 
the single atom qubit case [3, [ig as all operations must 
use blockade to prohibit multi-atom excitation. Gate op¬ 
erations are performed via the collective, singly excited 
Rydberg state 

1 ^ 

where \rj) is the Rydberg state of the atom. A single 
qubit rotation R{9, (jf) with area 9 and phase (f between 
ensemble states |0), |I) is implemented as the three pulse 
sequence |I) |f), |r) < > 10), If) 11). Note 

that the coupling strength between states |T), |f) is the 
single atom Rabi frequency U while the coupling between 
|0), |f) is at the collective Rabi frequency fljv = VN^- 
Since UIn depends on N, the one-qubit gate pulse lengths 
depend on the number of atoms. A Cz gate between 
control and target ensembles c, t is implemented as the 
three pulse sequence |I)c |f)c, |i)t -A- |f)t, |f)c ^ 

_ TT 27r TT 

|l)c. The Cz gate pulses do not depend on the number 
of atoms. The N dependence of the one-qubit gates can 
be strongly suppressed using adiabatic pulse sequences so 
that high fidelity gate operations are possible with small, 
but unknown values of A^[T^. 

The experimental setting is as described in [l8l |. In 
brief we prepare a cold sample of ®^Rb atoms in a 
magneto-optical trap (MOT) and then load a variable 
number of atoms into optical dipole traps. The dipole 
traps shown in Fig. [T] are formed by focusing 1064 nm 
light to waists (1/e^ intensity radii) of 3.0 p,m. The atoms 
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FIG. 1. (color online) Experimental geometry a) and transi¬ 
tions used for qubit control b). The Raman light is only used 
for preparation of product states, as discussed in connection 
with Fig. [3] 


are cooled to a temperature of ^ 150 /rK in 1-1.5 mK 
deep optical potentials. This gives approximately Gaus¬ 
sian shaped density distributions with typical standard 
deviations cr_L = 0.7 /rm perpendicular to the long trap 
axis and Uz = 1 /rm parallel to the long axis. The es¬ 
timated density at trap center is n/N = 5 x 10^® m“^. 
We apply a bias magnetic field along the trap axis of 
Bz = 0.24 mT and optically pump into |0) = |5si/2,/ = 
2,771/ = 0) using TT polarized 795 nm light resonant with 
|5 si/ 2, / = 2) —>• |5pi/2, / = 2) and 780 nm repump light 
resonant with |5si/2,/ = 1) -)■ |5p3/2,/ = 2). 

Rydberg excitation coupling |0),|f) is performed by 
off-resonant two-photon transitions via 5 pq /t[ 19| using 
counter-propagating 780o and 480 nm light. With (j+ 
polarization for both beams we couple to the Rydberg 
state |r) = \nd^/ 2 ,'rnj = 5/2) which is selected with a 
Bz = 0.37 mT bias field. The other qubit ground state 
is |1) = |5 si/ 2,/ = l,m/ = 0). Coupling between |I), |r) 
is performed with 780i and 480 nm light where 780o and 
780i have the same propagation vector and polarization 
but a frequency difference of 6.8 GHz corresponding to 
the ®^Rb / = 1 O / = 2 clock frequency. In the experi¬ 
ments reported below we used Rydberg levels 97^5/2 and 
In both cases strong blockade was observed in 
individual ensembles with no evidence for double excita¬ 
tion of the logical |1) state[l^. While we do not observe 
double excitation of |1), experiments with two ensembles 
do show evidence for double excitation of the Rydberg 
state |f), which plays a role in limiting the fidelity with 
which we can prepare the |1) state. 

We proceed to demonstrate the coherence of the en¬ 
semble states of Eq. © using Ramsey interferometry. 
The amplitude of the Ramsey signal is used to quantify 
the presence of N atom entanglement in the ensemble, 
as has been observed in other recent experiments [ill. . 
Details of the analysis showing that 82 ± 6% of the atoms 
participate in the entangle^W) state are presented in 
the supplemental material [22|. We load 3 < iV < 10 



FIG. 2. (color online) Ramsey interference measurement of 
qubit coherence for N = 7.6. The peak-peak amplitude of the 
oscillation as a function of the gap time gives = 2.6(3) ms. 
The circles are data points with Ter error bars and the dashed 
and solid lines are fits to the functions Ua(t),Ub(t) defined in 
the text. The gap time is the time t between the i?i(7r) pulses 
in Eq. ©• All data have been corrected for ~ 1.5% probabil¬ 
ity per atom of the blow away giving an unwanted transition 
from |0) ^ |1). The inset shows the Ramsey oscillations for 
gap times of 0 (solid line), 0.5 ms (dashed line), and 2.5 ms 
(dashed-dotted line). 

atoms into one of the optical traps. The number of atoms 
loaded for each measurement follows a Poisson distribu¬ 
tion with mean N. Each measurement starts with optical 
pumping into |0) followed by the pulse sequence 

It/) = Ri(7r)Ro(^/2)i?iWG(t)Ri(7r)Ro(^/2)|0). (2) 

Here Ro{9) is a pulse of area 9 between states |0), |r) and 
Ri{9) is a pulse of area 9 between states |1), |7). The 
first Rq{'k/2) pulse creates an equal superposition . 

This is then mapped to with a i?i(7r) pulse, we 

wait a gap time t described by an operator G(t), map 
|1) —>■ |f) with a i?i(7r) pulse, and then perform another 
tt/2 pulse between |0),|f). Finally, any population left 
in |r) is mapped back to |1) with another Ri(tt) pulse. 
Atoms in state |0) are then pushed out of the trap us¬ 
ing unbalanced radiation pressure from a beam resonant 
with |5 si/ 2,/ = 2) —^ |5p3/2i/ = 3) while the dipole 
trap light is chopped on and off. For the push out step a 
bias field is applied along x the narrow axis of the dipole 
traps, and the circularly polarized push out beam prop¬ 
agates along X. This is followed by a measurement of 
the number of atoms remaining in the dipole trap. The 
resulting data are shown in Fig. [5] The amplitude of the 
Ramsey interference at short gap times is limited by the 
I W) state preparation fidelity of about 50% for the atom 
number used in the figure. The fidelities of the Rq^tt) and 
Ri{tt) pulses used to prepare \W) are estimated to each 
be at least 90% on the basis of previous experiments [l^ 
and the strong inter-ensemble blockade effect we report 
below. We attribute the limited \W) state preparation 
fidelity to Rydberg dephasing, as will be discussed in 
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FIG. 3. (color online) Dependence of ensemble coherence time 
on N for \W) states (red circles) and product states (blue 
squares). The horizontal error bars represent the bounds 
for atom number measurements interleaved between Ramsey 
measurements. The open symbols are for preselected = 1 
states. The dashed lines are a guide to the eye. 
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the following. Periodic fluorescence measurements of the 
mean atom number (described in the supplemental ma¬ 
terial to [l^) bound drifts to 6.7 < TV < 9, during the 12 
hour measurement of this data set. 

The principal sources of decoherence in this experi¬ 
ment are expected to be ma gne tic noise, motional de¬ 
phasing, and atomic collisions [23j. For small atom num¬ 
bers and low collision rates we fit the Ramsey signal to 
the expression [ 2 ^ v\^{t,T 2 ) = fo/[l + 
and in the collision dominated regime we use a Gaus¬ 
sian form Va,{t) = where uq is the amplitude 

at t = 0. Both functional forms give the same T 2 time 
within our experimental error bars of T 2 = 2.6 ± 0.3 ms. 
The TT pulse times were 0.24 ^s for |0) —>■ |f), 0.06 ^s for 
the gap between pulses, and 0.68 ^s for |r) —>• | 1 ) giving a 
coherence to i?( 7 r) gate time ratio of approximately 2600. 

To further clarify the sensitivity to collisional dephas¬ 
ing Fig. [3] shows the measured T 2 for different N, in¬ 
cluding the case of = 1 Fock states which are selected 
using an additional fluorescence measurement before the 
Ramsey sequence [I^. We see that T 2 ~ UN in contrast 
to the 1/iV^ scaling observed for GHZ states [ 2 ^. The ob¬ 
served 1/7V scaling for \ W) states is expected for decoher¬ 
ence dominated by collisions since the collision rate per 
atom is proportional to N. For comparison, the T 2 time 
was also measured for product states |^/>) ^ (| 0 ) —i|l))®'^. 
These states were prepared using a two-frequency Raman 
laser coupling | 0 ) and | 1 ) via the 6 ^ 3/2 level [ 2 ^ as shown 
in Fig. [TJ Gomparison of the |1) {\W) state) and product 
state coherence data suggests that for iV > 5 the coher¬ 
ence time is limited by collisions. For TV < 5 as well as for 
the TV = 1 Fock state data the product states show longer 
coherence time. The coherence of the \W) states is mea¬ 
sured by comparison with a phase reference defined by 
the beatnote of the 780o and 780i Rydberg lasers which 
have a measured beatnote linewidth of 100 Hz FWHM. 
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FIG. 4. (color online) Ensemble to ensemble blockade for 
TVc = 9.9, iVt = 6.2. a) Probability of preparing |I)t without 
blockade (red circles, solid line) and with blockade (black cir¬ 
cles, dashed line). The solid line is a fit to a decaying sinusoid 
function from [1^. The dashed line is the same fit scaled by 
11%. b) Blockade data post selected on detection of |I)c. The 
dashed-dotted lines in both panels show the expected signal 
due to state leakage during blow-away in the control and tar¬ 
get regions. 


This linewidth is consistent with the observed shorter co¬ 
herence time of the \W) states compared to the product 
states which are referenced to the Raman laser beatnote 
which is in turn locked to a stable 6.8 GHz microwave 
oscillator. We anticipate that compensated optical traps 
and dynamical decoupling methods together with an op¬ 
tical lattice to reduce collisional effects can be used to 
greatly extend these coherence times [ 2 ^. 

To demonstrate ensemble-ensemble blockade we load 
atoms into control (c) and target (t) dipole traps, opti¬ 
cally pump into |0)c|0)t and apply one of two sequences. 
Preparation of a superposition of |0) and |1) in the 
target qubit is effected by the sequence [/a| 0 )c| 0 )t = 
i?i,t( 7 >')T?o,t( 0 )|O)c|O)t. This should ideally leave the 
qubits in the joint state |0)c [cos(0/2)|O)t — sin(d/2)|l)t] 
with the probability of preparing |I)t proportional to 
sin^(d/2), as is shown in Fig. 0^). We see the expected 
time dependence with a peak probability of ~ 

0.52, consistent with our earlier study of Fock state 
preparationfisj. 

Rydberg blockade between two ensembles 
is observed with the sequence C7b|0)c|0)t = 

i?i,c( 7 r)T?i,t( 7 r)T?o,t(d)i?o,c( 7 r)| 0 )c| 0 )t. Here we have 
used state | 0 ) of the control ensemble to block the target 
transfer with the final Ri.cCtt) pulse ideally leaving 
the qubits in the joint state |l)c|0)t. The data in Fig. 
gji) show a ratio of P\i),t{Uh)/P\i)^t{Ua) = 0 . 11 ( 1 ), 
i.e. a blockade fidelity of 0.89. This implies that the 
success probability of the transition T?o,c(7r)|0)c —t \ f)c is 
bounded below by the |l)t population ratio for the two 
sequences. We infer that at least one atom is excited to 
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the Rydberg state |r)c with probability > 0.89(1). 

As a further check on the inter-site blockade fidelity, 
events where the control site ends in state |I)c after se¬ 
quence Ub are post selected. The observed post-selected 
target population is shown in Figure |3 }d), along with the 
expected blow-away leakage rate of the control and target 
sites which is measured to be 0.2%/atom. From the data 
it can be seen that the post-selected results are consistent 
with perfect inter-site blockade. 

The observed high blockade fidelity exceeds that origi- 
nally achieved in experiments with single atom qubits 28, 
[ 29 II , and is certainly sufficient to create entanglement be¬ 
tween ensemble qubits. What has so far limited a demon¬ 
stration of deterministic entanglement is the relatively 
low probability of up to 62% [18| with which the ensem¬ 
ble state 11) can be prepared. In order to gain insight into 
what is limiting the state preparation fidelity we looked 
for signatures of Rydberg-Rydberg interactions concur¬ 
rently with strong blockade. Ideally the probability of 
preparing |I)c with sequence Ub, should be independent 
of the pulse area 9 applied to the target ensemble. How¬ 
ever a clear dependence on 6 can be seen in Fig. [S^). 
We believe this effect is due to long range interactions, 
where the amplitude for Rydberg atom excitation in the 
target site is sufficiently blockaded to prevent it from 
making the transfer to |l)t with any significant proba¬ 
bility, yet the target ensemble Rydberg excitation still 
interacts with the control ensemble strongly enough to 
disrupt the control ensemble state transfer. A similar 
situation of partial blockade together with decoherence 
of multi-atom ground-Rydberg Rabi oscillations was re¬ 
ported earlier in [l9| . 

A two-atom Rydberg interaction effect should scale 
with the Rydberg double excitation probability, i.e. P-?. oc 
where B is the ensemble mean blockade shift [30j. 
To check this, we extract the slopes from linear fits to the 
data for small 9 and compare to the scaling pa¬ 
rameter 


F — 

^ - “at. 


{R/RoY\ 


-2 


^ I^double- 


(3) 


Here n is the Rydberg principal quantum number and R 
is the site - site separation. The larger F is for a given 
set of parameters, the stronger the Rydberg-Rydberg in¬ 
teraction, and thus the larger the slope of 
Indeed, this is the behavior we observe, as shown in Fig. 
Eb), for a range of N, R, and n. 

This interaction effect hints at the possible mecha¬ 
nism responsible for the observed reduction in the prob¬ 
ability P|i) of preparing the collective qubit state in a 
single ensemble. The spatial extent of one ensemble is 
~ 2az = 14 ^m giving a length scale in between the 
lower two data sets in Fig. [5^). The intra-ensemble Ryd¬ 
berg interactions are significantly stronger than between 
atoms located in different ensembles at the same sepa¬ 
ration because the dipole-dipole interaction angular fac¬ 



FIG. 5. (color online) Probability of preparing state |l)c as a 
function of the target ensemble pulse area 6. a) Probability 
for several parameter sets: (Illd 5 y 2 , R = 8.3 and 8.7 fim) 
(red diamonds), (97^5/2, R = 8.3 and 8.7 pm) (green circles), 
(97d5/2, R = 17 pm) (yellow squares). The data has been 
normalized to 1 at 0 = 0 for clarity, with typical success 
probability 40-60%. b) Comparison of the slope of the data 
in panel (a) with the scaling parameter F from Eq. ([Sll. The 
color markers are the same as in panel a). 


tors favor atom pairs separated along z[30j. These con¬ 
siderations imply that lack of perfect blockade leading 
to long range Rydberg-Rydberg interactions in a single 
ensemble only partially explains the observed maximum 
of P|i) = 0.62 18|. Another candidate explanation is 
very strong interactions at short range in a single en¬ 
semble which mix levels together and open anti-blockade 
resonance channels 3l| . The doubly excited molecular 
energy structure becomes difficult to calculate with con¬ 
fidence at short range, with many molecular potentials 
near resonant [s^. For our typical Rydberg state 97^5/2 
this characteristic separation is ^ 5 pm, and for a 6 atom 
sample with our ensemble spatial distributions an average 
of 7 atom pairs out of 15 have R < 5 pm. We conjecture 
that the strong, short range interactions give an ampli¬ 
tude for double excitation, resulting in Rydberg-Rydberg 
interactions which dephase the ground-Rydberg rotations 
needed for state preparation, thereby limiting the prob¬ 
ability of preparing the ensemble |1) state. A related re¬ 
duction of the fidelity of Rydberg mediated atom-photon 
coupling in dense ensembles due to Rydberg-ground state 
interactions has also been observed [Uj. 

In conclusion, we have demonstrated the coherence of 
ensemble qubit basis states. The coherence time scales 
approximately inversely with the number of atoms, but 
is still several ms and 2600 times longer than our char¬ 
acteristic gate time for N ^ 10. Additionally we have 
demonstrated inter-ensemble blockade with a fidelity of 
0.89 and ^ 1.0 when post-selecting on control ensem¬ 
ble excitation. We identified Rydberg-Rydberg interac¬ 
tions from weak double excitations, either at long or short 
range, as a possible mechanism limiting the fidelity of en¬ 
semble state preparation. Future work towards ensemble 
entanglement and quantum computation will explore the 
use of a background optical lattice to better localize the 
ensembles while limiting uncontrolled short range inter¬ 
actions. 

This work was funded by NSF grant PHY-1104531 and 
the AFOSR Quantum Memories MURI. 
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Coherence and Rydberg blockade of atomic en¬ 
semble qubits 
Supplementary Material 


MULTIPARTITE W-STATE ENTANGLEMENT 
VERIFICATION 

In order to demonstrate multipartite entanglement it 
is necessary to show that the results obtained in a mea¬ 
surement cannot be reproduced with a separable state. 
Thus we require that the iV-particle state in question 
\ip^) satisfies 

(4) 

for any K in the range N/2 < K < N. In this sup¬ 
plemental material we evaluate the observed signatures 
of IE-state entanglement. These signatures include the 
•\/]V-enhancement of the Rabi frequency between |0) and 
|f), and the amplitude of the Ramsey oscillations. 

Collective Rabi Frequency Enhancement 

The interaction of an ensemble with a light field can 
be written in the basis of individual atom excitations 
|{0,1}^^^) 0 |{0,1}^^^) 0 1}^^^). The Hamiltonian 

Hint describing the evolution of the system is a block 
tridiagonal matrix. The basis states are denoted as \nk), 
where 0 < n < V is the eigenvalue of the excitation 
number operator M = J2k=i the index 

k labels the degenerate eigenstates, e.g. |Ii) = |10 • • • 0), 
II 2 ) = |0I • • • 0), etc. . Here is the effective 

spin operator for atom k along z. In this basis Hint is 
given by: 

Hint = A -|- A = 

Aq A(o4) 0 • • • 0 ■ 

'^{N-2,N-l) ^N-1 ■^{N-1,N) 

- 0 0 Apf 

(5) 

The matrix Hint has dimensions 2^ x 2^, the dimension 
of N 2-level systems. The dimension of the block di¬ 
agonal sub-matrices is given by the binomial coefficient, 
dim(A„) = (^) = Nn- The sub-matrices, A„, contain 
information concerning the sub-systems specific energy 
levels 

Nn 

A„ = (6) 


where <5^"^ refers to the energy of the fc**' basis state in 
the subspace with eigenvalue n. 

The matrices on the upper and lower diagonals couple 
states with excitation numbers differing by ±1, \nk) 

|n ± Ij) with coupling strength a defined by 

= ajk\{n±l)j){nk\ 

= \in±l):j){{n±l)j\A\nk){nkl (7) 

where 

N 

A=Y, (8) 

m—1 

and am is the strength of the light-atom coupling at atom 
m. In an ideal Rydberg blockaded ensemble states with 
n > 1 are not excited and all am are equal. Departures 
from the ideal case are accounted for by allowing for atom 
specific am and double excitations are included by trun¬ 
cating the basis at n = 2 and adding the doubly excited 
interaction energies to A 2 . 

A strong blockade shift, 6m~'^'^ = Sdd ^ ctm reduces 
the available Hilbert space for the problem to n = {0,1}, 
and Hint becomes: 


0 

Ol 

02 ■ 

■ OiN 


ai 


0 • 

■ 0 


02 

0 

r(l) 

0 

(9) 

Oat 

0 

0 • 

1 

■■ 



The detunings Sm'^ are nominally 0, so it makes sense 
to treat the entries as a perturbation. Under the 
condition of perfect blockade and no detuning, the en¬ 
ergy eigenstates of Hint = A are the dressed states 
(|0) ± 11)) with total angular momentum J = N/2 
and N — 1 orthogonal states with total angular mo¬ 
mentum J = {N/2 - I): (|0) ± |I)), |(I)_l)}, where 

|i) = Efci with a% = The eigen¬ 

values determine the speed at which the system evolves, 
for (|0) ± |1)) this speed is ia^v implying a collective 

enhancement of -s/Na when the coupling strengths are 
homogeneous. 

Our system has low inhomogeneous coupling contribu¬ 
tions as evidenced by the oln = scaling ob¬ 

served in our previous work [l8| , for reference an average 
scaling of 0.972 is predicted from experimental param¬ 
eters. The observation of a/A scaling of the coupling 
strength is a classic signature of Rydberg blockade and N 
participating wavefunctions, as the |1) state is the only 
state that can evolve with that coupling strength. A 
state with /c-partite entanglement consistent with the ob¬ 
served perfect blockade given by, \tp^) = |I^') 0 10-^“^), 
will still oscillate at the same Vn frequency, but the 
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FIG. 6. Monte Carlo calculations for N = 5 atoms for 
10,000 randomized instances of atom positions and velocities 
consistent with our experimental parameters of (a) Projection 
of a symmetric eigenstate without inhomogeneous broaden¬ 
ing, | —(|0} — |1)), onto the energy eigenstate of the 
full inhomogeneous Hamiltonian Hint I-} = E-\-) and (b) 
projection of the symmetric eigenstate onto the orthogonal 
subspace {1(1)*^')}. 

amplitude will be reduced to the overlap with 
|(1^|'0^)P = k/N, this is discussed further in the next 
section. 

The orthogonal singly-excited states |(T)_l) do not cou¬ 
ple to the symmetric states {|0), |I)} under ideal condi¬ 
tions (5^^^ = 0). This becomes clear when the Bloch 
picture is invoked, since the symmetric states have to¬ 
tal angular momentum J = N/2 while the |(i)_L) states 
have J = N/2 — 1 and a rotation on the Bloch sphere 
conserves angular momentum. Inhomogeneous broaden¬ 
ing, including differential AC Stark shifts, Doppler shifts, 
and finite intermediate state lifetimes, are added pertur- 
batively with A and provide a mechanism for coupling 
into the |(1 )_l) space. This coupling should be negligi¬ 
ble and reduce with increasing N and additionally will 
not display the characteristic VN enhancement. Figure 
inishows simulated projections of (|0) — |1)) along the 


energy eigenstates of Hint ior N = 5 atoms with our 
experimental parameters. 

Coherence Amplitude 

Since the coupling to the orthogonal subspace is neg¬ 
ligible for our experimental parameters, the amplitude 
of the Ramsey fringe oscillations provide a threshold for 
entanglement. A thermal sample of singly excited states 
|lth) = Sfci where c/k is a random phase fac¬ 

tor for the fc*** atom, will only couple back to | 0 ) by the 
amount of overlap with the |1) state. The projection 
|(i|lth)P will average to 1/A, therefore an oscillation 
with contrast above 1/A cannot be a thermal sample. 

To generate a threshold for /c-partite entanglement we 
perform a numerical simulation along the lines of the 
analysis in 0 - Briefly, the goal is to generate an upper 
bound on a measurement of Pi = |(I|'i/)p as a function of 
^0 = K 0 |V')P for states I'lp) with a maximum of k entan¬ 
gled particles. We establish bounds in two ways. First, 
we do not assume Rydberg blocakde so multiple excita¬ 
tions are possible. This is done by creating a random 
fc-partite entangled wavefunction 

IV') = -IV’m-l) (10) 

where = sin ( 0 ^/ 2 ) -b cos ( 0 ^/ 2 ) e*'^‘ 6i 

and (pi are randomly generated, and km = A — (m— l)k. 
We extract the maximum Pi for a given Pq bin obtained 
numerically to arrive at the thresholds shown in Fig. [Tli) 
for fc = 3 particle entanglement with ensemble atom num¬ 
bers A = 4 — 8. Any state above the threshold must 
have at least fc-partite entanglement. The black cross is 
an experimental data point recorded for a sample with 
A = 8.8 atoms, verifying the presence of entanglement. 

Rydberg blockade limits the Hilbert space to n < 1 ex¬ 
citations, which simplifies the calculation and enables an 
analytical bound for the fc-partite entanglement thresh¬ 
old. The state in Eq. (flUl) includes kets with multiple 
excitations. To remove these we impose the blockade 
condition P(„>i) = 0 and write the state as 

\tp) = (ailoW) + 0 (11) 

Maximization of Pj for a given Pg can be readily accom¬ 
plished analytically to give 

^i = ^(l-^o)- ( 12 ) 

Note that this a.grees with the limiting case of 
Pq = 0 from [20[ . Rearranging m to give 

U pmax 

, and given our extreme value 
(Pq,Pj) = (0.44 ± 0.02,0.46 ± 0.03) we can show 
that we meet the threshold for creation of the VF-state 































FIG. 7. (color online) a) Numerically determined bounds 
for {Po,Pi) using Eq. (IIUII . Rydberg blockade is not as¬ 
sumed in the calculation so multiple excitations are allowed. 
States above the {N, k) line imply there are at least k en¬ 
tangled atoms in the N atom ensemble. Calculated bounds 
for = 4 — 8 are shown, top to bottom. The dashed black 
line shows the amplitudes for the thermal singly-excited state 
11th) with N = A. The solid black line represents the range of 
experimental Ramsey oscillation data with the cross showing 
the value at tgap = 0 ms from Fig. 2 in the main text, b) An¬ 
alytical bounds assuming perfect blockade using Eqs. (Unnj. 
The entanglement thresholds are the straight lines shown for 
N = 9 and fc = 9 — 6, top to bottom. The data shown by the 
black line and cross exceeds the k = 7 threshold. 


with ^ = 82 ± 6%. Similar arguments for the presence 
of entanglement based on the amplitude of Ramsey 
oscillations have been used in 211. 


SUMMARY 


In summary we have shown evidence for N particle 
W-state entanglement on the basis of the following three 
arguments. First, the excellent agreement of the ob¬ 
served collectively enhanced Rabi frequency with theory 
reported in our previous work[l8j| using the same experi¬ 
mental setup and procedures as are used here implies an 
N component wavefunction. Second, the amplitude of 
the Ramsey-style oscillations for N = 8.8 is four times 
larger than the 1 /N limit expected from a thermal sam¬ 
ple of singly excited states. Our data shows entangle¬ 
ment without making the assumption of perfect block¬ 
ade. Third, with the assumptions of perfect blockade, 
entanglement percentage independent of iV, and negligi¬ 
ble coupling to {|(1)_l)}, which is justified by Fig. 1, then 
^ = 82 ±6%. In other words 82 ±6% of the atoms in the 
ensemble are participating in the FF-state entanglement. 
This result is not changed in a statistically significant way 
when compared with simulations based on experimental 
parameters that include imperfect blockade. 












